Rockfish Historical Harvest Estimates

This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.

Issue	Howard	Bayes
Time series	1999 - present	1977 - present
Bias in SWHS	Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests.	Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases	Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data.	Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations	Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas.	Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation	Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with.	By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Results

• Model: LB_fit_3pH

• Model Run: Gen4int_indcomp_swhsR_FULL_pHB4pars_re0full_altwt_slepssquz_thru2024_6e+05_2025-11-09 (Naming convention: modelname_thruYEAR_iters_rundate)

• Proportion of parmaters converged (Rhat <= 1.1): 0.9975

• See end of document for summary of unconverged parameters.

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

And black rockfish releases…

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Total Biomass Removal Estimates

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Model fit

Logbook residuals

Figure 20.- Residuals from logbook harvests.

SWHS residuals

Figure 21.- Residuals from SWHS harvests.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

Figure 23.- Mean percent of harvest by charter anglers.

When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(slope|non-pelagic & non-yellowye) For release estimates

## NULL

Weight Fits

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.5046296

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.15 while 0.85 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable	area_n	year	user	mean	median	lower	upper	Rhat	sd	skew	area
p_pelagic	11	1	1	0.467	0.467	0.228	0.712	1.011	0.126	0.000	CSEO
p_pelagic	11	11	2	0.404	0.399	0.186	0.652	1.015	0.121	0.041	CSEO
p_pelagic	11	15	1	0.430	0.426	0.198	0.684	1.017	0.124	0.032	CSEO
p_pelagic	11	15	2	0.394	0.389	0.172	0.646	1.015	0.122	0.041	CSEO
p_pelagic	11	17	2	0.401	0.395	0.182	0.640	1.014	0.120	0.050	CSEO
p_pelagic	11	2	2	0.400	0.397	0.183	0.652	1.013	0.121	0.025	CSEO
p_pelagic	11	22	2	0.410	0.403	0.192	0.659	1.013	0.120	0.058	CSEO
p_pelagic	11	4	1	0.468	0.466	0.232	0.709	1.012	0.126	0.016	CSEO
p_pelagic	11	5	2	0.399	0.392	0.183	0.656	1.011	0.123	0.057	CSEO
p_pelagic	11	7	1	0.468	0.466	0.236	0.713	1.012	0.125	0.016	CSEO
p_pelagic	11	8	1	0.467	0.466	0.228	0.712	1.011	0.127	0.008	CSEO
p_pelagic	11	9	1	0.467	0.467	0.230	0.713	1.017	0.125	0.000	CSEO
p_pelagic	13	1	1	0.301	0.281	0.101	0.595	1.014	0.132	0.152	NSEI
p_pelagic	13	1	2	0.123	0.107	0.035	0.313	1.014	0.071	0.224	NSEI
p_pelagic	13	10	2	0.131	0.115	0.037	0.308	1.016	0.073	0.220	NSEI
p_pelagic	13	11	2	0.131	0.114	0.039	0.320	1.019	0.074	0.230	NSEI
p_pelagic	13	13	2	0.137	0.121	0.041	0.315	1.013	0.073	0.219	NSEI
p_pelagic	13	14	1	0.310	0.298	0.119	0.566	1.010	0.119	0.101	NSEI
p_pelagic	13	15	1	0.284	0.268	0.106	0.554	1.010	0.117	0.137	NSEI
p_pelagic	13	15	2	0.131	0.117	0.041	0.295	1.010	0.067	0.208	NSEI
p_pelagic	13	2	1	0.300	0.281	0.098	0.585	1.015	0.131	0.145	NSEI
p_pelagic	13	2	2	0.124	0.105	0.033	0.312	1.014	0.073	0.259	NSEI
p_pelagic	13	3	1	0.301	0.282	0.100	0.598	1.015	0.132	0.144	NSEI
p_pelagic	13	4	1	0.302	0.283	0.108	0.588	1.010	0.128	0.148	NSEI
p_pelagic	13	4	2	0.125	0.110	0.033	0.310	1.018	0.073	0.204	NSEI
p_pelagic	13	5	2	0.124	0.107	0.034	0.311	1.012	0.072	0.236	NSEI
p_pelagic	13	6	2	0.126	0.109	0.034	0.315	1.012	0.073	0.231	NSEI
p_pelagic	13	7	2	0.126	0.111	0.035	0.308	1.015	0.072	0.210	NSEI
p_pelagic	13	8	1	0.308	0.292	0.108	0.598	1.012	0.128	0.125	NSEI
p_pelagic	13	9	2	0.129	0.111	0.036	0.308	1.015	0.073	0.248	NSEI
p_pelagic	16	1	1	0.400	0.392	0.158	0.673	1.032	0.138	0.058	SSEO
p_pelagic	16	1	2	0.376	0.369	0.142	0.654	1.033	0.136	0.052	SSEO
p_pelagic	16	10	1	0.403	0.397	0.160	0.676	1.040	0.137	0.044	SSEO
p_pelagic	16	10	2	0.382	0.372	0.149	0.657	1.033	0.135	0.074	SSEO
p_pelagic	16	11	1	0.404	0.400	0.163	0.685	1.041	0.136	0.029	SSEO
p_pelagic	16	11	2	0.385	0.379	0.147	0.660	1.034	0.134	0.045	SSEO
p_pelagic	16	12	1	0.401	0.397	0.159	0.670	1.028	0.135	0.030	SSEO
p_pelagic	16	12	2	0.384	0.377	0.152	0.661	1.022	0.133	0.053	SSEO
p_pelagic	16	13	1	0.405	0.398	0.162	0.674	1.030	0.135	0.052	SSEO
p_pelagic	16	13	2	0.383	0.376	0.148	0.646	1.024	0.132	0.053	SSEO
p_pelagic	16	14	1	0.391	0.384	0.156	0.655	1.027	0.130	0.054	SSEO
p_pelagic	16	14	2	0.382	0.373	0.160	0.648	1.030	0.130	0.069	SSEO
p_pelagic	16	15	1	0.396	0.391	0.163	0.657	1.034	0.131	0.038	SSEO
p_pelagic	16	15	2	0.386	0.376	0.156	0.656	1.027	0.133	0.075	SSEO
p_pelagic	16	16	1	0.391	0.382	0.163	0.659	1.022	0.130	0.069	SSEO
p_pelagic	16	16	2	0.384	0.377	0.155	0.654	1.017	0.130	0.054	SSEO
p_pelagic	16	17	1	0.446	0.441	0.208	0.710	1.020	0.131	0.038	SSEO
p_pelagic	16	17	2	0.399	0.392	0.171	0.663	1.037	0.131	0.053	SSEO
p_pelagic	16	18	1	0.426	0.419	0.192	0.684	1.020	0.129	0.054	SSEO
p_pelagic	16	18	2	0.399	0.395	0.170	0.663	1.022	0.128	0.031	SSEO
p_pelagic	16	19	1	0.420	0.414	0.188	0.687	1.015	0.130	0.046	SSEO
p_pelagic	16	19	2	0.399	0.394	0.170	0.659	1.022	0.128	0.039	SSEO
p_pelagic	16	2	1	0.400	0.394	0.148	0.682	1.028	0.138	0.043	SSEO
p_pelagic	16	2	2	0.378	0.372	0.137	0.644	1.035	0.135	0.045	SSEO
p_pelagic	16	20	1	0.414	0.410	0.190	0.665	1.019	0.124	0.032	SSEO
p_pelagic	16	20	2	0.404	0.398	0.175	0.671	1.013	0.128	0.047	SSEO
p_pelagic	16	21	1	0.398	0.389	0.180	0.659	1.014	0.124	0.073	SSEO
p_pelagic	16	21	2	0.399	0.393	0.184	0.655	1.016	0.125	0.048	SSEO
p_pelagic	16	22	2	0.405	0.399	0.187	0.664	1.014	0.124	0.048	SSEO
p_pelagic	16	3	1	0.398	0.395	0.149	0.676	1.037	0.139	0.022	SSEO
p_pelagic	16	3	2	0.376	0.369	0.137	0.650	1.037	0.136	0.051	SSEO
p_dsr	16	31	1	0.407	0.407	0.376	0.437	1.016	0.016	0.000	SSEO
p_pelagic	16	4	1	0.404	0.403	0.155	0.677	1.031	0.141	0.007	SSEO
p_pelagic	16	4	2	0.379	0.378	0.136	0.648	1.029	0.136	0.007	SSEO
p_pelagic	16	5	1	0.402	0.397	0.148	0.670	1.040	0.138	0.036	SSEO
p_pelagic	16	5	2	0.378	0.371	0.140	0.659	1.034	0.136	0.051	SSEO
p_pelagic	16	6	1	0.399	0.396	0.149	0.678	1.030	0.139	0.022	SSEO
p_pelagic	16	6	2	0.380	0.371	0.142	0.654	1.029	0.135	0.067	SSEO
p_pelagic	16	7	1	0.402	0.399	0.152	0.677	1.045	0.139	0.022	SSEO
p_pelagic	16	7	2	0.376	0.370	0.137	0.656	1.034	0.135	0.045	SSEO
p_pelagic	16	8	1	0.405	0.400	0.154	0.670	1.043	0.137	0.036	SSEO
p_pelagic	16	8	2	0.380	0.373	0.148	0.655	1.038	0.135	0.052	SSEO
p_pelagic	16	9	1	0.403	0.399	0.155	0.682	1.038	0.140	0.029	SSEO
p_pelagic	16	9	2	0.383	0.380	0.137	0.660	1.036	0.138	0.022	SSEO
p_pelagic	4	1	1	0.537	0.548	0.261	0.788	1.010	0.137	-0.080	PWSO
p_pelagic	4	1	2	0.544	0.553	0.271	0.795	1.012	0.138	-0.065	PWSO
p_pelagic	4	10	1	0.536	0.541	0.271	0.780	1.013	0.135	-0.037	PWSO
p_yellow	4	10	1	0.878	0.893	0.694	0.975	1.011	0.075	-0.199	PWSO
p_pelagic	4	11	1	0.536	0.541	0.266	0.780	1.011	0.136	-0.037	PWSO
p_yellow	4	11	2	0.875	0.890	0.693	0.975	1.011	0.073	-0.204	PWSO
p_pelagic	4	12	2	0.545	0.552	0.267	0.791	1.012	0.137	-0.051	PWSO
p_yellow	4	12	1	0.878	0.893	0.689	0.972	1.016	0.072	-0.208	PWSO
p_pelagic	4	14	1	0.531	0.536	0.267	0.772	1.012	0.131	-0.038	PWSO
p_pelagic	4	14	2	0.543	0.548	0.265	0.796	1.014	0.135	-0.037	PWSO
p_yellow	4	14	1	0.878	0.892	0.702	0.975	1.024	0.074	-0.190	PWSO
p_pelagic	4	15	2	0.545	0.549	0.269	0.782	1.014	0.132	-0.030	PWSO
p_yellow	4	15	1	0.878	0.896	0.690	0.971	1.012	0.073	-0.246	PWSO
p_pelagic	4	16	1	0.534	0.540	0.268	0.789	1.010	0.136	-0.044	PWSO
p_yellow	4	16	2	0.875	0.892	0.685	0.971	1.013	0.074	-0.230	PWSO
p_pelagic	4	17	2	0.549	0.556	0.279	0.782	1.010	0.133	-0.053	PWSO
p_pelagic	4	19	1	0.537	0.538	0.275	0.781	1.013	0.131	-0.008	PWSO
p_pelagic	4	2	1	0.535	0.540	0.249	0.781	1.016	0.136	-0.037	PWSO
p_pelagic	4	2	2	0.542	0.547	0.273	0.790	1.014	0.134	-0.037	PWSO
p_pelagic	4	3	2	0.542	0.548	0.260	0.792	1.011	0.138	-0.043	PWSO
p_yellow	4	3	1	0.880	0.895	0.696	0.979	1.012	0.073	-0.205	PWSO
p_pelagic	4	4	1	0.537	0.539	0.265	0.776	1.014	0.136	-0.015	PWSO
p_yellow	4	4	1	0.881	0.897	0.700	0.979	1.010	0.073	-0.218	PWSO
p_pelagic	4	5	1	0.535	0.541	0.261	0.786	1.011	0.135	-0.044	PWSO
p_pelagic	4	5	2	0.541	0.545	0.255	0.791	1.014	0.138	-0.029	PWSO
p_pelagic	4	6	1	0.536	0.541	0.261	0.784	1.013	0.137	-0.037	PWSO
p_pelagic	4	6	2	0.545	0.552	0.258	0.791	1.011	0.136	-0.051	PWSO
p_pelagic	4	7	1	0.537	0.540	0.265	0.786	1.018	0.137	-0.022	PWSO
p_pelagic	4	7	2	0.544	0.550	0.265	0.783	1.011	0.135	-0.045	PWSO
p_pelagic	4	8	1	0.540	0.542	0.265	0.792	1.010	0.137	-0.015	PWSO
p_yellow	4	8	2	0.876	0.890	0.686	0.975	1.011	0.075	-0.186	PWSO
p_pelagic	4	9	2	0.544	0.549	0.283	0.785	1.013	0.134	-0.037	PWSO
p_black	8	40	1	0.461	0.414	0.340	0.957	1.021	0.160	0.293	afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable	area_n	year/species	mean	median	lower	upper	Rhat	sd	skew	year	area
beta0_pH	11	4	0.783	0.914	-1.487	1.902	1.268	0.764	-0.171		CSEO
beta1_pH	11	4	2.010	1.918	0.715	4.221	1.213	0.792	0.116		CSEO
beta0_pelagic	11		-0.137	-0.111	-0.604	0.148	1.130	0.188	-0.138		CSEO
beta0_pH	12	4	-0.146	0.914	-17.076	2.226	1.472	4.209	-0.252		EWYKT
beta1_pH	12	4	4.233	3.160	1.802	20.983	1.469	4.216	0.254		EWYKT
beta3_pH	12	4	41.902	42.116	40.206	43.764	1.156	0.939	-0.228		EWYKT
beta0_pH	13	4	0.401	0.736	-3.051	1.586	1.374	1.116	-0.300		NSEI
beta1_pH	13	4	1.593	1.290	0.313	5.004	1.383	1.118	0.271		NSEI
beta0_pH	14	4	0.310	0.435	-1.491	1.276	1.375	0.698	-0.179		NSEO
beta1_pH	14	4	2.056	1.942	0.905	3.976	1.363	0.769	0.148		NSEO
beta0_yellow	14		-0.298	-0.236	-1.341	-0.010	1.268	0.350	-0.177		NSEO
beta0_pH	15	4	0.486	0.684	-1.237	1.532	1.364	0.732	-0.270		SSEI
beta1_pH	15	4	1.515	1.359	0.323	3.261	1.306	0.765	0.204		SSEI
beta3_pH	15	4	41.309	41.266	38.831	43.706	1.203	1.200	0.036		SSEI
beta0_pH	16	4	0.761	0.907	-1.699	2.152	1.344	0.963	-0.152		SSEO
beta1_pH	16	4	3.619	3.460	2.104	6.231	1.324	1.013	0.157		SSEO
beta0_pelagic	16		-0.452	-0.361	-1.340	0.071	1.133	0.384	-0.237		SSEO
mu_beta2_pH	2	1	1.449	0.766	0.073	7.081	1.131	2.083	0.328		NG
mu_beta0_pH	3	4	0.442	0.773	-3.752	1.657	1.433	1.356	-0.244		PWSI
mu_beta1_pH	3	4	2.274	2.031	0.642	5.233	1.230	1.312	0.185		PWSI
beta2_pH	6	1	1.312	0.468	0.108	8.151	1.122	2.231	0.378		SOKO2SAP
beta2_pH	8	1	1.108	0.390	0.114	7.392	1.116	1.841	0.390		afognak
beta3_pH	9	1	31.955	31.625	27.153	39.698	1.114	2.817	0.117		eastside

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic	non.BSAI.SOKO.Harv	BSAI.SOKO.Harv	non.BSAI.SOKO.Rel	BSAI.SOKO.Rel
n unconverged at Rhat = 1.01	307.00	193.00	944.00	332.00
Prop. of total	0.61	0.39	0.74	0.26
Prop. median less than 1	0.32	0.23	0.35	0.45
Prop. median less than 10	0.46	0.57	0.58	0.80
Prop. median less than 50	0.67	0.78	0.76	0.86
Median of the median values	15.18	7.21	4.29	1.51
Mean of the median values	234.96	87.28	157.68	90.52
Mean skew of posteriors	0.26	0.30	0.27	0.24

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic	non.BSAI.SOKO.Harv	BSAI.SOKO.Harv	non.BSAI.SOKO.Rel	BSAI.SOKO.Rel
n unconverged at Rhat = 1.1	3.00	4.00	43.00	40.00
Prop. of total	0.43	0.57	0.52	0.48
Prop. median less than 1	1.00	1.00	0.79	0.52
Prop. median less than 10	1.00	1.00	0.86	0.92
Prop. median less than 50	1.00	1.00	0.95	0.95
Median of the median values	0.00	0.54	0.01	0.51
Mean of the median values	0.01	0.53	10.50	10.71
Mean skew of posteriors	0.15	0.20	0.14	0.17

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable	Freq
beta0_pelagic	2
beta0_pH	6
beta0_yellow	1
beta1_pH	6
beta2_pH	2
beta3_pH	3
eps_pH	3
Ho_ay	1
Ho_ayu	4
Hs_ayu	1
Hy_ay	3
Hy_ayu	3
mu_beta0_pH	1
mu_beta1_pH	1
mu_beta2_pH	1
pH	2
R_ayg	1
R_ayu	1
Rb_ayg	1
Rb_ayu	1
Rd_ayg	1
Rd_ayu	1
re_pH	5
Ro_ay	15
Ro_ayg	16
Ro_ayu	28
Rp_ayg	1
Rp_ayu	1
Rs_ay	2
Rs_ayg	8
Rs_ayu	19
Ry_ay	12
Ry_ayg	3
Ry_ayu	24

Parameter estimates:

Summary Table of Parameter Estimates

Parameter	mean	sd	Lower_CI	Median	Upper_CI
mu_bc_H[1]	-0.016	0.069	-0.142	-0.021	0.127
mu_bc_H[2]	0.072	0.055	-0.029	0.071	0.188
mu_bc_H[3]	-0.267	0.067	-0.398	-0.268	-0.135
mu_bc_H[4]	-0.924	0.180	-1.287	-0.918	-0.584
mu_bc_H[5]	1.313	1.102	0.120	1.063	4.002
mu_bc_H[6]	-2.400	0.313	-3.000	-2.413	-1.762
mu_bc_H[7]	-0.457	0.103	-0.662	-0.453	-0.266
mu_bc_H[8]	0.272	0.346	-0.308	0.241	1.067
mu_bc_H[9]	-0.247	0.125	-0.496	-0.248	0.005
mu_bc_H[10]	0.058	0.066	-0.065	0.057	0.191
mu_bc_H[11]	-0.126	0.037	-0.200	-0.125	-0.054
mu_bc_H[12]	-0.250	0.107	-0.487	-0.243	-0.053
mu_bc_H[13]	-0.160	0.072	-0.297	-0.162	-0.016
mu_bc_H[14]	-0.305	0.089	-0.488	-0.303	-0.127
mu_bc_H[15]	-0.339	0.047	-0.435	-0.339	-0.246
mu_bc_H[16]	-0.141	0.366	-0.785	-0.156	0.643
mu_bc_R[1]	1.042	0.121	0.803	1.042	1.276
mu_bc_R[2]	1.010	0.080	0.842	1.014	1.156
mu_bc_R[3]	0.224	0.090	0.035	0.225	0.399
mu_bc_R[4]	0.678	0.198	0.270	0.686	1.040
mu_bc_R[5]	1.674	0.375	0.933	1.673	2.392
mu_bc_R[6]	-1.429	0.339	-2.109	-1.426	-0.737
mu_bc_R[7]	0.509	0.171	0.157	0.510	0.836
mu_bc_R[8]	0.651	0.145	0.385	0.647	0.943
mu_bc_R[9]	0.383	0.167	0.036	0.390	0.694
mu_bc_R[10]	0.897	0.104	0.694	0.899	1.102
mu_bc_R[11]	1.113	0.071	0.981	1.111	1.255
mu_bc_R[12]	0.873	0.162	0.540	0.875	1.176
mu_bc_R[13]	1.031	0.089	0.847	1.031	1.199
mu_bc_R[14]	0.861	0.123	0.620	0.862	1.095
mu_bc_R[15]	0.801	0.097	0.603	0.802	0.987
mu_bc_R[16]	1.353	0.117	1.123	1.356	1.574
tau_pH[1]	2.425	0.183	2.091	2.421	2.806
tau_pH[2]	1.606	0.154	1.323	1.600	1.923
tau_pH[3]	1.196	0.113	0.989	1.191	1.429
tau_pH[4]	1.017	0.141	0.771	1.009	1.318
tau_pH[5]	0.427	0.076	0.295	0.422	0.594
beta0_pH[1,1]	0.549	0.183	0.175	0.553	0.917
beta0_pH[2,1]	1.460	0.176	1.102	1.468	1.772
beta0_pH[3,1]	1.389	0.226	0.860	1.411	1.769
beta0_pH[4,1]	1.601	0.230	1.073	1.615	1.997
beta0_pH[5,1]	-0.776	0.414	-1.829	-0.728	-0.120
beta0_pH[6,1]	-0.523	0.495	-1.875	-0.377	0.090
beta0_pH[7,1]	0.757	0.011	0.736	0.757	0.778
beta0_pH[8,1]	-0.680	0.353	-1.636	-0.629	-0.120
beta0_pH[9,1]	-0.351	0.355	-1.122	-0.299	0.179
beta0_pH[10,1]	0.295	0.209	-0.128	0.310	0.683
beta0_pH[11,1]	0.069	0.311	-0.622	0.087	0.623
beta0_pH[12,1]	0.629	0.190	0.228	0.640	0.971
beta0_pH[13,1]	-0.654	0.151	-0.970	-0.650	-0.368
beta0_pH[14,1]	0.013	0.181	-0.350	0.016	0.361
beta0_pH[15,1]	0.066	0.193	-0.376	0.087	0.385
beta0_pH[16,1]	-0.879	0.221	-1.344	-0.865	-0.477
beta0_pH[1,2]	2.562	0.085	2.396	2.562	2.737
beta0_pH[2,2]	3.079	0.022	3.037	3.079	3.123
beta0_pH[3,2]	3.625	0.041	3.547	3.625	3.705
beta0_pH[4,2]	2.937	0.019	2.899	2.937	2.975
beta0_pH[5,2]	4.550	1.393	2.737	4.268	8.261
beta0_pH[6,2]	3.327	0.152	3.037	3.323	3.633
beta0_pH[7,2]	2.033	0.049	1.939	2.033	2.131
beta0_pH[8,2]	2.823	0.064	2.698	2.823	2.951
beta0_pH[9,2]	3.952	0.206	3.559	3.946	4.364
beta0_pH[10,2]	3.996	0.131	3.753	3.990	4.265
beta0_pH[11,2]	1.885	0.009	1.868	1.885	1.903
beta0_pH[12,2]	1.794	0.069	1.663	1.792	1.932
beta0_pH[13,2]	2.496	0.026	2.445	2.495	2.548
beta0_pH[14,2]	1.803	0.024	1.757	1.803	1.849
beta0_pH[15,2]	2.545	0.017	2.513	2.545	2.577
beta0_pH[16,2]	2.766	0.025	2.716	2.766	2.815
beta0_pH[1,3]	1.096	0.409	0.150	1.161	1.694
beta0_pH[2,3]	2.126	0.018	2.091	2.126	2.161
beta0_pH[3,3]	2.692	0.024	2.646	2.692	2.739
beta0_pH[4,3]	2.980	0.035	2.911	2.979	3.051
beta0_pH[5,3]	3.480	2.328	0.939	2.901	9.887
beta0_pH[6,3]	1.366	0.376	0.714	1.336	2.168
beta0_pH[7,3]	0.860	0.050	0.763	0.860	0.954
beta0_pH[8,3]	0.255	0.295	-0.277	0.295	0.609
beta0_pH[9,3]	-2.278	1.059	-4.852	-2.132	-0.572
beta0_pH[10,3]	-0.518	0.873	-2.580	-0.395	0.843
beta0_pH[11,3]	1.128	0.507	0.089	1.159	2.038
beta0_pH[12,3]	1.003	0.661	-0.367	0.996	2.312
beta0_pH[13,3]	1.025	0.482	0.002	1.035	1.911
beta0_pH[14,3]	0.603	0.426	-0.333	0.642	1.337
beta0_pH[15,3]	0.729	0.404	-0.098	0.736	1.466
beta0_pH[16,3]	1.215	0.626	0.035	1.194	2.477
beta0_pH[11,4]	0.783	0.764	-1.487	0.914	1.902
beta0_pH[12,4]	-0.146	4.209	-17.076	0.914	2.226
beta0_pH[13,4]	0.401	1.116	-3.051	0.736	1.586
beta0_pH[14,4]	0.310	0.698	-1.491	0.435	1.276
beta0_pH[15,4]	0.486	0.732	-1.237	0.684	1.532
beta0_pH[16,4]	0.761	0.963	-1.699	0.907	2.152
beta0_pH[11,5]	2.971	0.149	2.699	2.962	3.281
beta0_pH[12,5]	-0.914	0.508	-1.921	-0.897	0.053
beta0_pH[13,5]	4.892	0.417	4.122	4.879	5.764
beta0_pH[14,5]	3.909	0.363	3.257	3.887	4.691
beta0_pH[15,5]	2.286	0.124	2.059	2.280	2.550
beta0_pH[16,5]	0.247	0.064	0.160	0.235	0.406
beta1_pH[1,1]	2.585	0.486	1.658	2.562	3.665
beta1_pH[2,1]	1.552	0.377	0.981	1.521	2.312
beta1_pH[3,1]	1.828	0.455	1.113	1.762	3.063
beta1_pH[4,1]	1.874	0.440	1.180	1.813	2.976
beta1_pH[5,1]	2.298	0.546	1.450	2.227	3.670
beta1_pH[6,1]	2.145	0.856	1.061	1.886	4.443
beta1_pH[7,1]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,1]	2.327	0.933	1.167	2.073	4.667
beta1_pH[9,1]	2.061	0.670	1.171	1.929	4.019
beta1_pH[10,1]	2.149	0.468	1.460	2.067	3.300
beta1_pH[11,1]	3.372	0.548	2.517	3.303	4.733
beta1_pH[12,1]	2.897	0.279	2.363	2.887	3.466
beta1_pH[13,1]	3.452	0.245	2.988	3.449	3.945
beta1_pH[14,1]	2.840	0.338	2.202	2.826	3.572
beta1_pH[15,1]	2.397	0.308	1.900	2.369	3.013
beta1_pH[16,1]	4.037	0.305	3.462	4.024	4.639
beta1_pH[1,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[2,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[3,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[4,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[5,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[6,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[7,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[9,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[10,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[11,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[12,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[13,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[14,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[15,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[16,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[1,3]	1.902	1.250	0.041	1.755	5.208
beta1_pH[2,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[3,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[4,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[5,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[6,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[7,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,3]	4.659	2.401	0.920	4.160	10.655
beta1_pH[9,3]	7.584	2.660	3.586	7.144	13.822
beta1_pH[10,3]	6.617	2.726	2.823	6.048	12.825
beta1_pH[11,3]	2.157	0.551	1.082	2.160	3.263
beta1_pH[12,3]	3.217	0.681	1.947	3.186	4.644
beta1_pH[13,3]	1.953	0.519	0.951	1.949	3.023
beta1_pH[14,3]	2.168	0.486	1.300	2.128	3.225
beta1_pH[15,3]	2.096	0.474	1.195	2.084	3.070
beta1_pH[16,3]	3.660	0.685	2.216	3.680	4.916
beta1_pH[11,4]	2.010	0.792	0.715	1.918	4.221
beta1_pH[12,4]	4.233	4.216	1.802	3.160	20.983
beta1_pH[13,4]	1.593	1.118	0.313	1.290	5.004
beta1_pH[14,4]	2.056	0.769	0.905	1.942	3.976
beta1_pH[15,4]	1.515	0.765	0.323	1.359	3.261
beta1_pH[16,4]	3.619	1.013	2.104	3.460	6.231
beta1_pH[11,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[12,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[13,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[14,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[15,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[16,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[1,1]	1.237	0.850	0.329	0.914	3.388
beta2_pH[2,1]	1.681	0.960	0.262	1.520	3.421
beta2_pH[3,1]	1.359	0.847	0.193	1.215	2.950
beta2_pH[4,1]	1.168	0.802	0.223	0.916	2.932
beta2_pH[5,1]	2.402	2.413	0.257	1.713	8.942
beta2_pH[6,1]	1.312	2.231	0.108	0.468	8.151
beta2_pH[7,1]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,1]	1.108	1.841	0.114	0.390	7.392
beta2_pH[9,1]	1.289	1.889	0.146	0.560	7.832
beta2_pH[10,1]	1.504	1.805	0.227	0.854	7.287
beta2_pH[11,1]	0.415	0.212	0.179	0.366	0.951
beta2_pH[12,1]	1.569	0.816	0.625	1.309	3.674
beta2_pH[13,1]	1.013	0.439	0.552	0.905	2.364
beta2_pH[14,1]	1.160	0.686	0.431	0.949	3.245
beta2_pH[15,1]	1.049	0.600	0.357	0.894	2.796
beta2_pH[16,1]	1.055	0.398	0.609	0.968	1.997
beta2_pH[1,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[2,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[3,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[4,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[5,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[6,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[7,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[9,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[10,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[11,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[12,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[13,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[14,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[15,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[16,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[1,3]	6.602	6.597	0.091	5.072	21.568
beta2_pH[2,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[3,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[4,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[5,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[6,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[7,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,3]	1.557	1.462	0.168	1.075	6.155
beta2_pH[9,3]	0.301	0.605	0.087	0.200	0.988
beta2_pH[10,3]	0.437	0.752	0.101	0.267	1.967
beta2_pH[11,3]	-6.262	4.915	-18.723	-4.918	-0.607
beta2_pH[12,3]	-7.095	4.843	-20.189	-5.877	-1.204
beta2_pH[13,3]	-6.532	4.865	-19.717	-5.244	-1.111
beta2_pH[14,3]	-6.664	4.868	-19.994	-5.205	-1.256
beta2_pH[15,3]	-6.770	4.816	-19.696	-5.453	-1.150
beta2_pH[16,3]	-6.594	4.704	-19.623	-5.110	-1.569
beta2_pH[11,4]	-5.831	4.673	-17.967	-4.368	-0.705
beta2_pH[12,4]	-6.030	4.558	-17.637	-4.641	-1.047
beta2_pH[13,4]	-6.010	4.502	-17.306	-4.694	-0.891
beta2_pH[14,4]	-5.892	4.435	-17.086	-4.560	-0.975
beta2_pH[15,4]	-5.892	4.638	-17.482	-4.416	-0.785
beta2_pH[16,4]	-6.069	4.384	-17.450	-4.666	-1.260
beta2_pH[11,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[12,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[13,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[14,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[15,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[16,5]	0.000	0.000	0.000	0.000	0.000
beta3_pH[1,1]	34.964	1.166	32.822	34.932	37.571
beta3_pH[2,1]	32.921	1.410	30.446	32.839	35.931
beta3_pH[3,1]	33.053	1.619	29.984	32.946	37.087
beta3_pH[4,1]	32.790	1.944	29.509	32.737	36.917
beta3_pH[5,1]	27.625	1.151	25.931	27.493	30.731
beta3_pH[6,1]	33.278	3.203	26.471	33.142	40.548
beta3_pH[7,1]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,1]	34.309	3.526	28.742	33.897	41.704
beta3_pH[9,1]	31.955	2.817	27.153	31.625	39.698
beta3_pH[10,1]	32.752	1.507	30.207	32.648	35.872
beta3_pH[11,1]	31.412	1.234	29.292	31.368	33.954
beta3_pH[12,1]	30.783	0.471	29.792	30.796	31.636
beta3_pH[13,1]	31.424	0.527	30.434	31.398	32.594
beta3_pH[14,1]	32.500	0.753	31.236	32.445	34.144
beta3_pH[15,1]	31.028	0.716	29.534	31.071	32.348
beta3_pH[16,1]	30.122	0.375	29.421	30.117	30.854
beta3_pH[1,2]	34.505	5.500	25.488	34.340	43.496
beta3_pH[2,2]	34.337	5.461	25.395	34.284	43.494
beta3_pH[3,2]	30.610	7.308	18.521	30.649	42.396
beta3_pH[4,2]	30.438	7.191	18.575	30.352	42.386
beta3_pH[5,2]	30.483	7.276	18.554	30.700	42.317
beta3_pH[6,2]	30.526	7.255	18.609	30.672	42.400
beta3_pH[7,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[9,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[10,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[11,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[12,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[13,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[14,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[15,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[16,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[1,3]	38.067	3.494	29.280	39.471	42.696
beta3_pH[2,3]	34.586	5.458	25.509	34.742	43.530
beta3_pH[3,3]	30.436	7.245	18.629	30.300	42.398
beta3_pH[4,3]	30.460	7.207	18.617	30.630	42.465
beta3_pH[5,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[6,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[7,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,3]	42.013	1.295	38.783	42.058	43.857
beta3_pH[9,3]	36.138	1.849	34.531	35.410	41.612
beta3_pH[10,3]	37.949	2.580	34.668	37.220	43.369
beta3_pH[11,3]	40.736	1.085	38.895	40.748	42.714
beta3_pH[12,3]	42.009	0.740	40.296	42.268	42.851
beta3_pH[13,3]	41.413	0.731	39.959	41.419	42.794
beta3_pH[14,3]	40.477	0.656	39.132	40.497	41.756
beta3_pH[15,3]	40.625	0.637	39.191	40.659	41.816
beta3_pH[16,3]	41.838	0.388	41.021	41.897	42.590
beta3_pH[11,4]	40.624	1.240	38.192	40.520	43.369
beta3_pH[12,4]	41.902	0.939	40.206	42.116	43.764
beta3_pH[13,4]	42.012	1.198	39.582	42.048	43.890
beta3_pH[14,4]	40.528	0.812	38.948	40.556	42.090
beta3_pH[15,4]	41.309	1.200	38.831	41.266	43.706
beta3_pH[16,4]	41.772	0.483	40.923	41.802	42.823
beta3_pH[11,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[12,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[13,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[14,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[15,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[16,5]	25.000	0.000	25.000	25.000	25.000
beta4_pH[1,1]	0.263	0.021	0.224	0.263	0.305
beta4_pH[2,1]	0.211	0.005	0.201	0.211	0.221
beta4_pH[3,1]	0.173	0.007	0.159	0.172	0.187
beta4_pH[4,1]	0.294	0.011	0.273	0.294	0.316
beta4_pH[5,1]	0.424	0.259	0.089	0.397	0.938
beta4_pH[6,1]	0.794	0.506	0.313	0.732	1.518
beta4_pH[7,1]	0.859	0.331	0.355	0.810	1.628
beta4_pH[8,1]	0.158	0.069	0.044	0.152	0.313
beta4_pH[9,1]	0.213	0.085	0.060	0.209	0.387
beta4_pH[10,1]	0.001	0.001	0.000	0.001	0.003
beta4_pH[11,1]	0.446	0.012	0.424	0.446	0.471
beta4_pH[12,1]	0.641	0.023	0.596	0.641	0.687
beta4_pH[13,1]	0.632	0.014	0.605	0.631	0.658
beta4_pH[14,1]	0.669	0.032	0.610	0.668	0.734
beta4_pH[15,1]	0.797	0.021	0.756	0.797	0.838
beta4_pH[16,1]	1.220	0.042	1.136	1.218	1.306
beta4_pH[1,2]	0.307	0.096	0.142	0.301	0.509
beta4_pH[2,2]	0.513	0.023	0.469	0.513	0.557
beta4_pH[3,2]	0.446	0.012	0.423	0.446	0.470
beta4_pH[4,2]	0.654	0.037	0.582	0.653	0.730
beta4_pH[5,2]	0.716	0.742	0.105	0.531	2.311
beta4_pH[6,2]	0.899	0.815	0.190	0.728	2.604
beta4_pH[7,2]	0.890	0.429	0.351	0.809	1.874
beta4_pH[8,2]	0.844	0.203	0.485	0.835	1.262
beta4_pH[9,2]	1.213	1.294	0.300	0.898	4.249
beta4_pH[10,2]	0.455	0.052	0.358	0.454	0.561
beta4_pH[11,2]	1.270	0.050	1.177	1.269	1.370
beta4_pH[12,2]	0.579	0.108	0.386	0.573	0.810
beta4_pH[13,2]	0.947	0.044	0.861	0.946	1.031
beta4_pH[14,2]	0.933	0.063	0.813	0.933	1.063
beta4_pH[15,2]	1.011	0.026	0.965	1.010	1.065
beta4_pH[16,2]	1.408	0.075	1.267	1.405	1.558
beta4_pH[1,3]	0.285	0.089	0.125	0.281	0.476
beta4_pH[2,3]	0.052	0.017	0.019	0.051	0.087
beta4_pH[3,3]	0.052	0.009	0.034	0.052	0.071
beta4_pH[4,3]	0.173	0.032	0.112	0.173	0.236
beta4_pH[5,3]	0.754	0.793	0.120	0.570	2.619
beta4_pH[6,3]	0.885	0.913	0.206	0.706	2.562
beta4_pH[7,3]	0.888	0.403	0.343	0.814	1.918
beta4_pH[8,3]	0.415	0.199	0.105	0.387	0.887
beta4_pH[9,3]	0.974	1.061	0.164	0.722	3.108
beta4_pH[10,3]	0.061	0.042	0.003	0.055	0.157
beta4_pH[11,3]	0.363	0.016	0.333	0.363	0.395
beta4_pH[12,3]	0.439	0.042	0.363	0.438	0.525
beta4_pH[13,3]	0.838	0.037	0.766	0.839	0.911
beta4_pH[14,3]	0.827	0.068	0.700	0.825	0.963
beta4_pH[15,3]	0.635	0.017	0.602	0.636	0.668
beta4_pH[16,3]	0.639	0.035	0.574	0.638	0.712
beta4_pH[1,4]	0.633	0.328	0.210	0.570	1.413
beta4_pH[2,4]	0.798	0.852	0.116	0.554	3.031
beta4_pH[3,4]	0.825	0.785	0.116	0.596	2.957
beta4_pH[4,4]	0.721	0.459	0.207	0.615	1.894
beta4_pH[5,4]	0.755	0.829	0.115	0.557	2.553
beta4_pH[6,4]	0.919	0.793	0.222	0.731	2.891
beta4_pH[7,4]	0.893	0.469	0.344	0.818	1.877
beta4_pH[8,4]	0.739	0.603	0.158	0.583	2.277
beta4_pH[9,4]	0.780	0.875	0.110	0.567	2.862
beta4_pH[10,4]	1.519	3.206	0.007	0.387	11.651
beta4_pH[11,4]	0.444	0.017	0.411	0.444	0.477
beta4_pH[12,4]	0.562	0.055	0.457	0.561	0.674
beta4_pH[13,4]	0.926	0.036	0.856	0.927	0.998
beta4_pH[14,4]	0.923	0.076	0.781	0.919	1.083
beta4_pH[15,4]	0.714	0.020	0.675	0.714	0.754
beta4_pH[16,4]	0.687	0.042	0.610	0.687	0.771
beta4_pH[1,5]	0.642	0.344	0.212	0.577	1.508
beta4_pH[2,5]	0.785	0.785	0.109	0.578	2.666
beta4_pH[3,5]	0.848	0.895	0.121	0.598	3.112
beta4_pH[4,5]	0.719	0.458	0.207	0.616	1.798
beta4_pH[5,5]	0.761	0.798	0.110	0.565	2.579
beta4_pH[6,5]	0.888	0.768	0.207	0.720	2.704
beta4_pH[7,5]	0.890	0.462	0.342	0.815	1.825
beta4_pH[8,5]	0.728	0.643	0.146	0.597	2.096
beta4_pH[9,5]	0.813	1.062	0.130	0.569	2.891
beta4_pH[10,5]	1.493	3.120	0.006	0.383	11.405
beta4_pH[11,5]	0.603	0.052	0.504	0.602	0.710
beta4_pH[12,5]	0.587	0.169	0.313	0.569	0.982
beta4_pH[13,5]	0.775	0.062	0.662	0.772	0.904
beta4_pH[14,5]	0.843	0.117	0.630	0.838	1.090
beta4_pH[15,5]	1.145	0.073	1.011	1.143	1.291
beta4_pH[16,5]	18.912	3.926	10.908	19.246	24.721
beta0_pelagic[1]	2.268	0.036	2.202	2.267	2.342
beta0_pelagic[2]	1.537	0.007	1.524	1.538	1.551
beta0_pelagic[3]	0.258	0.206	-0.256	0.296	0.571
beta0_pelagic[4]	0.175	0.287	-0.577	0.216	0.648
beta0_pelagic[5]	1.443	0.155	1.157	1.438	1.752
beta0_pelagic[6]	1.353	0.059	1.240	1.353	1.468
beta0_pelagic[7]	1.479	0.039	1.405	1.478	1.559
beta0_pelagic[8]	1.672	0.043	1.590	1.672	1.757
beta0_pelagic[9]	2.796	0.065	2.673	2.794	2.934
beta0_pelagic[10]	2.603	0.031	2.545	2.603	2.665
beta0_pelagic[11]	-0.137	0.188	-0.604	-0.111	0.148
beta0_pelagic[12]	1.801	0.036	1.734	1.800	1.875
beta0_pelagic[13]	-1.002	0.447	-1.706	-1.040	-0.095
beta0_pelagic[14]	-0.274	0.165	-0.664	-0.259	0.014
beta0_pelagic[15]	-0.210	0.040	-0.323	-0.195	-0.180
beta0_pelagic[16]	-0.452	0.384	-1.340	-0.361	0.071
beta1_pelagic[1]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[2]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[3]	0.690	0.230	0.299	0.671	1.268
beta1_pelagic[4]	0.878	0.318	0.355	0.845	1.716
beta1_pelagic[5]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[6]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[7]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[8]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[9]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[10]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[11]	1.140	0.248	0.731	1.111	1.764
beta1_pelagic[12]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[13]	2.782	0.607	1.452	2.794	3.880
beta1_pelagic[14]	1.781	0.430	1.123	1.707	2.834
beta1_pelagic[15]	0.077	0.133	0.000	0.013	0.472
beta1_pelagic[16]	1.531	0.554	0.743	1.433	2.798
beta2_pelagic[1]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[2]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[3]	7.045	6.279	0.280	5.413	22.823
beta2_pelagic[4]	7.268	6.306	0.232	5.863	22.633
beta2_pelagic[5]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[6]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[7]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[8]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[9]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[10]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[11]	1.848	1.496	0.308	1.350	5.705
beta2_pelagic[12]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[13]	0.185	0.079	0.095	0.165	0.393
beta2_pelagic[14]	0.868	0.863	0.208	0.596	3.513
beta2_pelagic[15]	2.628	1.914	0.156	2.216	7.115
beta2_pelagic[16]	0.817	1.025	0.131	0.401	4.188
beta3_pelagic[1]	28.785	3.903	22.343	28.826	35.207
beta3_pelagic[2]	28.828	3.855	22.444	28.838	35.186
beta3_pelagic[3]	28.353	2.572	23.630	29.087	32.278
beta3_pelagic[4]	25.942	1.511	22.626	26.198	28.957
beta3_pelagic[5]	29.459	6.670	18.630	29.406	40.487
beta3_pelagic[6]	29.512	6.639	18.517	29.434	40.372
beta3_pelagic[7]	29.462	6.604	18.511	29.222	40.516
beta3_pelagic[8]	29.512	6.758	18.453	29.568	40.566
beta3_pelagic[9]	29.460	6.798	18.527	29.520	40.545
beta3_pelagic[10]	29.618	6.716	18.498	29.707	40.537
beta3_pelagic[11]	30.107	1.192	27.755	30.045	32.697
beta3_pelagic[12]	30.000	0.000	30.000	30.000	30.000
beta3_pelagic[13]	31.011	4.080	24.808	30.302	40.218
beta3_pelagic[14]	34.285	1.717	31.420	34.121	38.255
beta3_pelagic[15]	33.842	5.745	23.479	34.289	41.666
beta3_pelagic[16]	31.052	2.461	26.257	31.096	35.934
beta0_yellow[1]	-0.485	0.079	-0.644	-0.485	-0.325
beta0_yellow[2]	0.541	0.082	0.367	0.550	0.678
beta0_yellow[3]	-0.212	0.129	-0.469	-0.217	0.014
beta0_yellow[4]	0.889	0.228	0.304	0.937	1.188
beta0_yellow[5]	0.253	0.328	-0.454	0.275	0.826
beta0_yellow[6]	0.815	0.108	0.601	0.815	1.020
beta0_yellow[7]	0.994	0.098	0.790	0.994	1.183
beta0_yellow[8]	1.026	0.094	0.843	1.027	1.210
beta0_yellow[9]	0.459	0.102	0.260	0.460	0.653
beta0_yellow[10]	0.540	0.071	0.399	0.540	0.675
beta0_yellow[11]	-0.204	0.129	-0.488	-0.194	0.020
beta0_yellow[12]	-1.081	0.093	-1.261	-1.081	-0.898
beta0_yellow[13]	-1.356	0.206	-1.779	-1.333	-1.056
beta0_yellow[14]	-0.298	0.350	-1.341	-0.236	-0.010
beta0_yellow[15]	-1.253	0.164	-1.603	-1.230	-1.032
beta0_yellow[16]	-0.243	0.093	-0.423	-0.244	-0.060
beta1_yellow[1]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[2]	1.014	0.234	0.582	1.013	1.488
beta1_yellow[3]	0.477	0.260	0.021	0.491	0.932
beta1_yellow[4]	1.329	0.603	0.643	1.174	3.155
beta1_yellow[5]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[6]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[7]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[8]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[9]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[10]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[11]	0.296	0.358	0.001	0.125	1.153
beta1_yellow[12]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[13]	1.808	1.666	0.045	1.065	6.006
beta1_yellow[14]	0.762	1.030	0.013	0.456	4.153
beta1_yellow[15]	0.170	0.232	0.011	0.076	0.796
beta1_yellow[16]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[1]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[2]	-7.546	5.711	-21.353	-6.328	-0.398
beta2_yellow[3]	-7.708	5.818	-22.184	-6.515	-0.429
beta2_yellow[4]	-5.507	5.803	-20.018	-3.484	-0.127
beta2_yellow[5]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[6]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[7]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[8]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[9]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[10]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[11]	-7.798	5.892	-22.129	-6.676	-0.361
beta2_yellow[12]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[13]	-5.331	5.231	-18.520	-3.379	-0.275
beta2_yellow[14]	-6.958	5.557	-20.018	-5.785	-0.121
beta2_yellow[15]	-7.451	5.438	-20.028	-6.408	-0.298
beta2_yellow[16]	0.000	0.000	0.000	0.000	0.000
beta3_yellow[1]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[2]	28.884	1.121	27.340	28.669	31.919
beta3_yellow[3]	32.388	3.304	23.384	32.568	38.907
beta3_yellow[4]	30.050	3.338	23.357	29.204	35.667
beta3_yellow[5]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[6]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[7]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[8]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[9]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[10]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[11]	35.065	3.671	30.314	34.078	43.333
beta3_yellow[12]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[13]	32.678	2.913	29.100	32.484	39.076
beta3_yellow[14]	33.971	3.617	29.237	33.465	42.953
beta3_yellow[15]	36.572	3.823	30.302	36.317	43.559
beta3_yellow[16]	35.000	0.000	35.000	35.000	35.000
beta0_black[1]	0.013	0.075	-0.133	0.012	0.160
beta0_black[2]	1.856	0.033	1.789	1.855	1.921
beta0_black[3]	1.237	0.037	1.166	1.237	1.311
beta0_black[4]	2.426	0.048	2.334	2.426	2.523
beta0_black[5]	5.188	2.352	1.906	4.716	11.142
beta0_black[6]	5.221	2.139	2.615	4.698	10.820
beta0_black[7]	4.095	2.188	1.612	3.509	10.131
beta0_black[8]	0.874	0.214	0.471	0.868	1.325
beta0_black[9]	2.802	0.269	2.283	2.798	3.345
beta0_black[10]	1.286	0.052	1.183	1.287	1.388
beta0_black[11]	2.979	0.051	2.878	2.979	3.078
beta0_black[12]	4.910	0.189	4.541	4.907	5.280
beta0_black[13]	0.436	0.333	0.018	0.353	1.223
beta0_black[14]	2.853	0.105	2.654	2.853	3.057
beta0_black[15]	0.723	0.053	0.620	0.723	0.825
beta0_black[16]	3.774	0.073	3.635	3.772	3.917
beta2_black[1]	0.000	0.000	0.000	0.000	0.000
beta2_black[2]	0.000	0.000	0.000	0.000	0.000
beta2_black[3]	0.000	0.000	0.000	0.000	0.000
beta2_black[4]	0.000	0.000	0.000	0.000	0.000
beta2_black[5]	0.000	0.000	0.000	0.000	0.000
beta2_black[6]	0.000	0.000	0.000	0.000	0.000
beta2_black[7]	0.000	0.000	0.000	0.000	0.000
beta2_black[8]	0.000	0.000	0.000	0.000	0.000
beta2_black[9]	0.000	0.000	0.000	0.000	0.000
beta2_black[10]	0.000	0.000	0.000	0.000	0.000
beta2_black[11]	0.000	0.000	0.000	0.000	0.000
beta2_black[12]	0.000	0.000	0.000	0.000	0.000
beta2_black[13]	-2.502	1.759	-6.994	-2.104	-0.286
beta2_black[14]	0.000	0.000	0.000	0.000	0.000
beta2_black[15]	0.000	0.000	0.000	0.000	0.000
beta2_black[16]	0.000	0.000	0.000	0.000	0.000
beta3_black[1]	25.000	0.000	25.000	25.000	25.000
beta3_black[2]	25.000	0.000	25.000	25.000	25.000
beta3_black[3]	25.000	0.000	25.000	25.000	25.000
beta3_black[4]	25.000	0.000	25.000	25.000	25.000
beta3_black[5]	25.000	0.000	25.000	25.000	25.000
beta3_black[6]	25.000	0.000	25.000	25.000	25.000
beta3_black[7]	25.000	0.000	25.000	25.000	25.000
beta3_black[8]	25.000	0.000	25.000	25.000	25.000
beta3_black[9]	25.000	0.000	25.000	25.000	25.000
beta3_black[10]	25.000	0.000	25.000	25.000	25.000
beta3_black[11]	25.000	0.000	25.000	25.000	25.000
beta3_black[12]	25.000	0.000	25.000	25.000	25.000
beta3_black[13]	37.402	3.242	27.764	38.103	40.941
beta3_black[14]	25.000	0.000	25.000	25.000	25.000
beta3_black[15]	25.000	0.000	25.000	25.000	25.000
beta3_black[16]	25.000	0.000	25.000	25.000	25.000
beta4_black[1]	-0.328	0.091	-0.502	-0.327	-0.152
beta4_black[2]	0.258	0.054	0.152	0.258	0.364
beta4_black[3]	-0.720	0.063	-0.843	-0.720	-0.594
beta4_black[4]	0.628	0.129	0.380	0.624	0.883
beta4_black[5]	0.794	1.502	-1.152	0.503	4.439
beta4_black[6]	0.897	1.814	-0.831	0.525	4.463
beta4_black[7]	0.776	1.448	-1.058	0.497	4.249
beta4_black[8]	0.086	0.279	-0.467	0.099	0.618
beta4_black[9]	1.036	0.810	-0.088	0.878	3.084
beta4_black[10]	0.073	0.068	-0.061	0.074	0.203
beta4_black[11]	-0.609	0.064	-0.729	-0.609	-0.479
beta4_black[12]	1.234	0.297	0.675	1.228	1.843
beta4_black[13]	-0.999	0.077	-1.156	-0.997	-0.854
beta4_black[14]	0.220	0.183	-0.128	0.218	0.586
beta4_black[15]	0.002	0.046	-0.089	0.002	0.094
beta4_black[16]	-0.418	0.116	-0.641	-0.419	-0.190
beta0_dsr[11]	-0.329	0.009	-0.348	-0.329	-0.312
beta0_dsr[12]	1.409	0.044	1.323	1.409	1.493
beta0_dsr[13]	0.155	0.021	0.115	0.155	0.197
beta0_dsr[14]	-0.127	0.017	-0.160	-0.128	-0.094
beta0_dsr[15]	0.205	0.014	0.178	0.206	0.232
beta0_dsr[16]	0.157	0.015	0.126	0.157	0.187
beta1_dsr[11]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[12]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[13]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[14]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[15]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[16]	0.000	0.000	0.000	0.000	0.000
beta4_dsr[11]	0.545	0.035	0.475	0.544	0.613
beta4_dsr[12]	0.802	0.232	0.353	0.802	1.261
beta4_dsr[13]	-0.079	0.039	-0.157	-0.079	-0.002
beta4_dsr[14]	-0.208	0.070	-0.347	-0.209	-0.072
beta4_dsr[15]	0.216	0.022	0.175	0.216	0.258
beta4_dsr[16]	-0.047	0.046	-0.138	-0.046	0.043
beta0_slope[11]	-2.208	0.021	-2.250	-2.208	-2.167
beta0_slope[12]	-4.989	0.186	-5.371	-4.982	-4.645
beta0_slope[13]	-1.501	0.033	-1.566	-1.500	-1.434
beta0_slope[14]	-2.712	0.051	-2.817	-2.711	-2.616
beta0_slope[15]	-1.592	0.022	-1.634	-1.592	-1.548
beta0_slope[16]	-2.649	0.043	-2.734	-2.651	-2.562
beta1_slope[11]	0.000	0.000	0.000	0.000	0.000
beta1_slope[12]	0.000	0.000	0.000	0.000	0.000
beta1_slope[13]	0.000	0.000	0.000	0.000	0.000
beta1_slope[14]	0.000	0.000	0.000	0.000	0.000
beta1_slope[15]	0.000	0.000	0.000	0.000	0.000
beta1_slope[16]	0.000	0.000	0.000	0.000	0.000
beta4_slope[11]	-0.336	0.075	-0.481	-0.337	-0.192
beta4_slope[12]	-4.641	0.776	-6.302	-4.597	-3.223
beta4_slope[13]	0.327	0.053	0.221	0.326	0.429
beta4_slope[14]	-0.160	0.152	-0.465	-0.157	0.127
beta4_slope[15]	-0.361	0.033	-0.424	-0.360	-0.295
beta4_slope[16]	-0.142	0.102	-0.350	-0.142	0.052
beta5_slope[11]	3.067	0.107	2.862	3.065	3.288
beta5_slope[12]	1.999	0.383	1.310	1.980	2.806
beta5_slope[13]	0.883	0.080	0.727	0.884	1.037
beta5_slope[14]	3.759	0.353	3.114	3.749	4.483
beta5_slope[15]	1.306	0.073	1.164	1.304	1.455
beta5_slope[16]	2.735	0.236	2.295	2.728	3.191
tau_comp[1]	3.766	0.300	3.201	3.754	4.403
tau_comp[2]	2.863	0.257	2.379	2.857	3.374
tau_comp[3]	0.870	0.074	0.733	0.868	1.016
tau_comp[4]	1.778	0.276	1.295	1.752	2.404
tau_comp[5]	1.262	0.189	0.917	1.253	1.665
sigma_H[1]	0.206	0.050	0.118	0.202	0.314
sigma_H[2]	0.175	0.028	0.127	0.172	0.237
sigma_H[3]	0.192	0.041	0.119	0.189	0.278
sigma_H[4]	0.461	0.080	0.328	0.452	0.644
sigma_H[5]	0.960	0.203	0.593	0.949	1.378
sigma_H[6]	0.520	0.187	0.160	0.511	0.909
sigma_H[7]	0.311	0.062	0.215	0.302	0.457
sigma_H[8]	0.417	0.087	0.274	0.408	0.608
sigma_H[9]	0.513	0.119	0.331	0.498	0.787
sigma_H[10]	0.211	0.039	0.144	0.206	0.298
sigma_H[11]	0.266	0.044	0.194	0.261	0.361
sigma_H[12]	0.382	0.147	0.203	0.333	0.724
sigma_H[13]	0.217	0.037	0.155	0.214	0.299
sigma_H[14]	0.499	0.090	0.346	0.490	0.698
sigma_H[15]	0.232	0.037	0.170	0.228	0.314
sigma_H[16]	0.216	0.041	0.150	0.212	0.309
lambda_H[1]	2.232	3.442	0.078	1.035	12.502
lambda_H[2]	5.018	4.694	0.460	3.615	17.385
lambda_H[3]	5.710	8.714	0.258	2.791	29.423
lambda_H[4]	0.006	0.004	0.001	0.005	0.018
lambda_H[5]	3.147	9.218	0.030	0.693	21.190
lambda_H[6]	6.287	14.043	0.009	0.207	46.884
lambda_H[7]	0.012	0.008	0.002	0.010	0.034
lambda_H[8]	7.949	10.190	0.118	4.414	36.763
lambda_H[9]	0.015	0.010	0.003	0.013	0.040
lambda_H[10]	0.228	0.363	0.026	0.156	0.839
lambda_H[11]	0.242	0.422	0.010	0.114	1.147
lambda_H[12]	3.674	4.879	0.167	2.096	18.229
lambda_H[13]	3.131	3.061	0.278	2.219	10.836
lambda_H[14]	3.476	4.156	0.249	2.185	14.993
lambda_H[15]	0.019	0.025	0.003	0.014	0.061
lambda_H[16]	0.466	0.748	0.031	0.246	2.331
mu_lambda_H[1]	4.068	1.902	1.000	3.871	8.133
mu_lambda_H[2]	3.674	1.957	0.507	3.467	7.813
mu_lambda_H[3]	3.382	1.852	0.683	3.116	7.600
sigma_lambda_H[1]	8.079	4.341	1.666	7.238	18.319
sigma_lambda_H[2]	8.138	4.737	0.794	7.512	18.273
sigma_lambda_H[3]	6.158	4.020	0.886	5.253	16.120
beta_H[1,1]	6.491	1.330	3.249	6.715	8.408
beta_H[2,1]	9.756	0.625	8.432	9.792	10.908
beta_H[3,1]	7.840	0.823	5.887	7.935	9.180
beta_H[4,1]	10.632	7.579	-4.885	10.613	25.119
beta_H[5,1]	-0.254	2.558	-5.644	-0.081	4.079
beta_H[6,1]	2.999	4.301	-7.438	4.430	8.226
beta_H[7,1]	1.467	5.674	-10.511	1.875	11.549
beta_H[8,1]	1.058	2.117	-2.661	1.185	3.379
beta_H[9,1]	13.674	5.582	2.746	13.617	24.629
beta_H[10,1]	6.847	1.854	2.967	6.913	10.294
beta_H[11,1]	4.992	3.542	-2.901	5.741	10.003
beta_H[12,1]	2.581	1.085	0.647	2.497	4.892
beta_H[13,1]	9.096	0.897	7.214	9.153	10.653
beta_H[14,1]	2.143	1.008	0.257	2.118	4.143
beta_H[15,1]	-7.051	3.380	-13.149	-7.154	0.155
beta_H[16,1]	3.983	2.979	-0.802	3.594	10.649
beta_H[1,2]	7.776	0.268	7.220	7.792	8.278
beta_H[2,2]	9.886	0.163	9.560	9.884	10.206
beta_H[3,2]	8.795	0.201	8.407	8.793	9.194
beta_H[4,2]	2.994	1.444	0.294	2.994	5.915
beta_H[5,2]	1.772	1.003	-0.145	1.764	3.739
beta_H[6,2]	6.039	1.107	3.536	6.191	7.800
beta_H[7,2]	2.377	1.073	0.494	2.306	4.623
beta_H[8,2]	3.053	0.784	1.474	3.137	4.198
beta_H[9,2]	3.153	1.104	1.000	3.133	5.399
beta_H[10,2]	8.028	0.369	7.267	8.039	8.719
beta_H[11,2]	9.814	0.638	8.857	9.705	11.273
beta_H[12,2]	3.938	0.373	3.221	3.927	4.699
beta_H[13,2]	9.162	0.245	8.695	9.155	9.653
beta_H[14,2]	4.031	0.343	3.393	4.018	4.749
beta_H[15,2]	11.556	0.612	10.287	11.565	12.663
beta_H[16,2]	4.463	0.834	2.924	4.418	6.112
beta_H[1,3]	8.444	0.264	7.965	8.435	8.982
beta_H[2,3]	9.865	0.150	9.576	9.863	10.172
beta_H[3,3]	9.455	0.167	9.145	9.451	9.793
beta_H[4,3]	-2.167	0.849	-3.897	-2.172	-0.572
beta_H[5,3]	3.900	0.629	2.602	3.912	5.126
beta_H[6,3]	8.406	1.148	6.745	8.151	10.760
beta_H[7,3]	-2.566	0.627	-3.870	-2.559	-1.372
beta_H[8,3]	5.266	0.419	4.675	5.203	6.307
beta_H[9,3]	-2.698	0.725	-4.140	-2.709	-1.296
beta_H[10,3]	8.553	0.277	8.022	8.549	9.108
beta_H[11,3]	8.503	0.295	7.868	8.527	9.010
beta_H[12,3]	5.265	0.321	4.512	5.306	5.809
beta_H[13,3]	8.858	0.176	8.496	8.859	9.199
beta_H[14,3]	5.781	0.275	5.181	5.803	6.255
beta_H[15,3]	10.283	0.298	9.702	10.281	10.879
beta_H[16,3]	5.877	0.595	4.668	5.914	6.943
beta_H[1,4]	8.101	0.185	7.714	8.110	8.442
beta_H[2,4]	9.850	0.134	9.574	9.855	10.104
beta_H[3,4]	9.947	0.154	9.605	9.961	10.223
beta_H[4,4]	11.815	0.465	10.881	11.826	12.714
beta_H[5,4]	5.585	0.746	4.357	5.500	7.311
beta_H[6,4]	7.166	0.917	5.284	7.359	8.478
beta_H[7,4]	8.551	0.344	7.884	8.551	9.227
beta_H[8,4]	6.726	0.227	6.264	6.732	7.140
beta_H[9,4]	7.429	0.449	6.560	7.418	8.336
beta_H[10,4]	7.526	0.224	7.099	7.518	7.982
beta_H[11,4]	9.511	0.192	9.123	9.516	9.888
beta_H[12,4]	7.248	0.192	6.881	7.249	7.623
beta_H[13,4]	9.067	0.140	8.784	9.071	9.338
beta_H[14,4]	7.834	0.218	7.428	7.831	8.272
beta_H[15,4]	9.473	0.221	9.034	9.473	9.899
beta_H[16,4]	9.435	0.239	8.977	9.429	9.912
beta_H[1,5]	8.894	0.158	8.572	8.901	9.186
beta_H[2,5]	10.679	0.101	10.493	10.675	10.896
beta_H[3,5]	10.755	0.164	10.456	10.748	11.086
beta_H[4,5]	8.347	0.490	7.411	8.326	9.335
beta_H[5,5]	5.189	0.597	3.699	5.262	6.151
beta_H[6,5]	8.981	0.614	8.051	8.870	10.308
beta_H[7,5]	6.577	0.341	5.910	6.569	7.273
beta_H[8,5]	8.241	0.193	7.917	8.223	8.643
beta_H[9,5]	8.187	0.462	7.249	8.186	9.111
beta_H[10,5]	10.047	0.230	9.579	10.048	10.490
beta_H[11,5]	11.594	0.221	11.162	11.591	12.036
beta_H[12,5]	8.610	0.190	8.240	8.606	8.990
beta_H[13,5]	10.142	0.130	9.888	10.138	10.414
beta_H[14,5]	9.234	0.216	8.820	9.227	9.695
beta_H[15,5]	11.315	0.233	10.858	11.315	11.769
beta_H[16,5]	9.828	0.186	9.437	9.835	10.162
beta_H[1,6]	10.123	0.244	9.756	10.085	10.715
beta_H[2,6]	11.468	0.117	11.235	11.468	11.692
beta_H[3,6]	10.620	0.155	10.287	10.629	10.899
beta_H[4,6]	12.923	0.867	11.191	12.936	14.588
beta_H[5,6]	5.669	0.627	4.625	5.622	6.982
beta_H[6,6]	8.720	0.576	7.419	8.752	9.781
beta_H[7,6]	10.566	0.589	9.351	10.570	11.702
beta_H[8,6]	9.534	0.246	9.061	9.539	9.967
beta_H[9,6]	8.289	0.768	6.795	8.263	9.850
beta_H[10,6]	9.136	0.340	8.429	9.159	9.762
beta_H[11,6]	10.688	0.359	9.917	10.720	11.308
beta_H[12,6]	9.368	0.241	8.887	9.360	9.871
beta_H[13,6]	11.094	0.154	10.801	11.091	11.412
beta_H[14,6]	9.921	0.286	9.308	9.932	10.453
beta_H[15,6]	10.550	0.417	9.713	10.548	11.344
beta_H[16,6]	10.583	0.260	10.040	10.594	11.066
beta_H[1,7]	10.134	1.265	6.780	10.423	11.883
beta_H[2,7]	12.123	0.521	10.980	12.145	13.090
beta_H[3,7]	10.299	0.671	8.783	10.355	11.394
beta_H[4,7]	1.401	4.503	-7.282	1.372	10.710
beta_H[5,7]	5.850	2.049	2.148	5.899	10.031
beta_H[6,7]	7.052	3.178	-0.977	8.150	10.735
beta_H[7,7]	7.597	2.995	1.751	7.562	13.794
beta_H[8,7]	10.854	0.838	9.293	10.861	12.350
beta_H[9,7]	4.959	3.930	-3.126	5.097	12.647
beta_H[10,7]	10.057	1.568	7.252	9.934	13.456
beta_H[11,7]	11.345	1.747	8.314	11.190	15.288
beta_H[12,7]	9.782	0.954	7.724	9.870	11.381
beta_H[13,7]	11.676	0.708	10.086	11.761	12.828
beta_H[14,7]	10.633	0.896	8.807	10.674	12.271
beta_H[15,7]	14.004	2.201	9.819	14.025	18.470
beta_H[16,7]	12.911	1.403	10.484	12.806	16.084
beta0_H[1]	8.948	18.927	-28.629	8.875	50.028
beta0_H[2]	10.679	8.115	-5.825	10.708	27.227
beta0_H[3]	9.873	10.128	-10.216	9.958	30.575
beta0_H[4]	10.122	180.032	-353.269	9.747	382.268
beta0_H[5]	3.774	27.192	-47.567	4.353	50.618
beta0_H[6]	7.946	51.207	-103.158	7.583	125.900
beta0_H[7]	3.632	135.887	-266.884	2.863	275.821
beta0_H[8]	6.446	21.694	-16.589	6.445	27.981
beta0_H[9]	11.208	124.940	-237.541	9.693	263.576
beta0_H[10]	8.316	38.449	-73.297	8.837	83.369
beta0_H[11]	8.201	50.439	-99.031	8.776	119.188
beta0_H[12]	6.624	12.174	-17.587	6.492	31.183
beta0_H[13]	10.427	11.635	-11.569	10.116	32.925
beta0_H[14]	7.204	11.170	-15.766	7.250	30.174
beta0_H[15]	6.138	112.883	-218.945	4.213	239.389
beta0_H[16]	7.333	33.649	-62.087	7.973	71.790